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Abstract 

We investigated the back reaction of cosmological perturbations on the evo- 
lution of the universe using the second order perturbation of the Einstein's 
equation. To incorporate the back reaction effect due to the inhomogeneity 
into the framework of the cosmological perturbation, we used the renormal- 
ization group method. The second order zero mode solution which appears 
by the non-linearities of the Einstein's equation is regarded as a secular term 
of the perturbative expansion, we renormalized a constant of integration con- 
tained in the background solution and absorbed the secular term to this con- 
stant. For a dust dominated universe, using the second order gauge invariant 
quantity, we derived the renormalization group equation which determines 
the effective dynamics of the Friedman-Robertson- Walker universe with the 
back reaction effect in a gauge invariant manner. We obtained the solution 
of the renormalization group equation and found that perturbations of the 
scalar mode and the long wavelength tensor mode works as positive spatial 
curvature, and the short wavelength tensor mode as radiation fluid. 
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I. INTRODUCTION 



Our universe seems to be very close to a Friedman-Robertson- Walker (FRW) spacetime 
at length scale of the order of the Hubble radius, but the metric and matter content appears 
to be highly inhomogeneous at smaller scale. The conventional cosmological perturbation 
approach |l|J^] treat such a situation as the homogeneous isotropic background plus small 
perturbation, and investigate the evolution of linear fluctuations. We must go beyond linear 
approximation to treat non-linear structure and construct suitable model of inhomogeneous 
universe which is close to a FRW universe on large scales. 

Based on the perturbative approach, the metric of inhomogeneous universe is approxi- 
mated by perturbative series and expanded as the background FRW metric, the first order 
linear perturbation and the second order perturbation which carries the non-linearities of the 
Einstein's equation. The non-linear effect affects not only the evolution of the fluctuation, 
but also change the evolution of the background FRW universe. This is the cosmological 
back reaction effect. Non- linear structure of the universe can affect the dynamics of the 
background FRW universe and it is important task to understand the effect of the cosmo- 
logical back reaction theoretically to determine the correct cosmological parameter of our 
universe. 

There are many works treating the back reaction problem so far. Issacson || considered 
the gravitational wave propagating on a slowly varying background spacetime. Assum- 
ing high frequency of the gravitational wave, he expanded the vacuum Einstein's equation 
perturbatively. The leading order of the expansion gives the propagation equation of the 
gravitational wave and the next order gives the equation which determines the slowly vary- 
ing background spacetime on which the gravitational wave propagates. The leading and the 
next order equations must be solved self-consistent way. Futamase [|,^f discussed the effect 
of the inhomogeneity on the global expansion rate of the universe using the cosmological 
post-Newtonian approximation with specific gauge conditions. Russ et al. @ treated the 
problem using the second order perturbation with synchronous comoving gauge and derived 
the evolution equation for the comoving volume. Boersma |7j used the linear perturbation 
with a uniform Hubble slice and derived the spatially averaged Einstein's equation. Abramo 



JJ derived the gauge invariant effective energy momentum tensor for the scalar field and 
the gravitational wave, and discussed the effect of inhomogeneity on the background FRW 
universe in a gauge invariant manner. 

These works are based on perturbation expansion of the Einstein's equation with the 
averaging of the metric. The averaging operation of the Einstein's equation is introduced 
to obtain the evolution equation of the effective scale factor. Papers f^-|7| use the different 
gauge conditions and the different perturbation scheme, and the relation of their result is 
not obvious. Abramo |8l-|il"H constructed the gauge invariant effective energy momentum 



tensor and discussed the effect of inhomogeneity on the background FRW universe, but he 
did not derived explicit form of the evolution equation of the effective scale factor. 

In this paper, we investigate the cosmological back reaction problem using the second 



order cosmological perturbation and the renormalization group method |[2|-[l9|]. We treat 
a dust dominated universe and examine how the inhomogeneity affect the evolution of the 
background FRW universe. We does not use the averaging procedure for the Einstein's 
equation to obtain the evolution equation of the effective scale factor. We only use the solu- 
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tion of the second order perturbation. Within the framework of the standard cosmological 
perturbation, we do not know how to include the second order non-linear effect in the back- 
ground FRW metric because the perturbation of each order separate. We use the method of 
the renormalization to incorporate the back reaction effect in the perturbation expansion. 
Using this method, we can interpret the back reaction as renormalization of constants of 
integration contained in the background FRW solution. Furthermore, it is easy to obtain 
the approximate analytic solution of the effective scale factor which contains the back reac- 
tion effect of inhomogeneity. To circumvent the problem of gauge ambiguity, we construct 
the second order gauge invariant quantities and evaluate the effect of the back reaction in 
a gauge invariant manner. We find that the scalar mode and the long wavelength tensor 
mode perturbation work as positive spatial curvature and the short wavelength tensor mode 
perturbation works as radiation fluid. 

The plan of this paper is as follows. In Sec. II, we introduce the renormalization group 
method by using FRW model. In Sec. Ill, we consider the second order cosmological 
perturbation of the Einstein's equation and obtain the zero mode solution of the second 
order perturbation. In Sec. IV, the second order gauge invariant quantity is introduced. 
The renormalization group method is applied to the solution of the second order zero mode 
perturbation and the effect of the back reaction is discussed. Sec. V is devoted to summary 
and discussion. We use the units in which c = h = 8nG = 1 throughout the paper. 



II. RENORMALIZATION GROUP METHOD 



In this section, we introduce the renormalization group method |Til2"|-iT)| using a FRW 
cosmological model. To obtain temporal evolution of the solution of the non-linear differ- 
ential equation, we usually apply a perturbative expansion. But naive perturbation often 
yields secular terms due to resonance phenomena. The secular term prevents us from getting 
approximate but global solutions. The renormalization group method is one of techniques 
to circumvent the problem. Starting from a naive perturbative expansion, the secular di- 
vergence is absorbed into constants of integration contained in the zeroth order solution by 
the renormalization procedure. The renormalized constants obey the renormalization group 
equation. 

We consider a dust dominated FRW universe with perfect fluid. The Einstein's equations 

are 

•2 1 1 

a = 3P0 + g Pi, (la) 
3 1 

a + -a 2 = --p!, (lb) 

where a is logarithm of a scale factor of the universe, po is the energy density of the dust 
field, pi and pi are the energy density and the pressure of perfect fluid with the equation of 
the state pi = (r — l)pi where T is a constant. Conservation equation for fluids are 

po + 3 a po = 0, (2a) 
pi + 3rdp!=0, (2b) 
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and their solutions are 



Po 



3 c 



Pi 



3 c\ 



where c ,Ci are constants of integration. Substitute this solution to Eq. ([Ta]) , we obtain 



a 



C Ci 



; >3o 



3ra 



(3) 



(4) 



We solve this equation perturbatively by assuming the second term of the right hand side 
is small: 

a = «o + oci + • • • . (5) 
The solution of the scale factor up to the first order of perturbation is given by 

3 c\ 



a(t) = a t 2/3 



1 + 



4 (3 - 2r) 



3I7 + 2-2r _ -f-2-2T\ 



(t 



(6) 



where ao and to are constants of integration of the zeroth order and the first order, respec- 
tively. For < r < 1, the first order term grows in time and the perturbation breaks 
down for large t. We improve this behavior by applying the renormalization group method 



We redefine the zeroth order integration constant a as 

a = a* (/■*) +Sa (t ;fx), 



(?) 



where \i is a renormalization point and Sao is a counter term which absorbs the secular 
divergence of the solution. The naive solution (Q) can be written 



a{t) = t 2 ' 3 

= t 2/3 



q„ ( n R\l-3T 

a«(n) + 5a«(t ;n) + 6Cl[ao) 



4 (3 - 2r) 



yU 2-2r + /i 2- 2 r_ t 2- 2 r ) 



%W+ 4(3 _ 2r) I* 



2-2V\ 



(8) 



where we have chosen the counter term 5a so as to absorb the (/i — / 
term: 



2_2r - t 2 ~ 2r ) dependent 



Sa (t ;ii) = o£(t ) - a* (aO 



3 Cl (a*) 1 -* 1 



4 (3 - 2f) 



(A* 



2-2r +2-2r\ 
l J 



This defines the renormalization transformation 
Kn,t ■ a§ ; (t ) 1 — > Oq (//) = a* (*o) - 



3 Cl (a*) 1 " 31, 
4 (3 - 2r) 



(// 



2-2r j.2-2r 







(9) 



(10) 



and this transformation forms the Lie group up to the first order of the perturbation. So we 

,2-2r tp~ 2r ) 

Lie group, and this makes it possible to produce a globally uniform approximated solution 



can have ao(/f) for arbitrary large value of (/z 2 21 — i$ ZL ) by assuming the property of the 
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of the original equation. The renormalization group equation is obtained by differentiating 
Eq.(||) with respect to fi, and setting t = /i: 







3ci(l T) i_ 3 r i-2r 

U/Q fjj 



dfi uu ~ y 2(3-2r) 
The renormalized solution is obtained by equating fi = t in 

9ciT 



(11) 



c + 



4 (3 - 2T) 



2-2F 



1/(3T) 



(12) 



Now we compare the renormalized solution with the naive solution for T = 2/3 and 
T = 4/3 cases. For T = 2/3, the energy density of the fluid is p\ oc a~ 2 and the effect of the 
fluid is equivalent to spatial curvature. The naive and the renormalized solutions are 



a(t) = a t 2/3 



1 + 



9ci 
20 a 2 



a R (t)=t^[c+ 9 -^t 2 / 3 



(t 2/3 -tT) 

1/2 



(13a) 
(13b) 



To see how the renormalized solution improves the naive solution, we choose = 1, C\ = 
— l,c = l,to = 0. Ci < corresponds to the positive spatial curvature. In this case, the 
exact solution of Eq. (01) is given by 
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COS Tj), 



(14) 



t = — (j) — sin rj) (0 < rj < 2ty). 



The scale factor of the exact solution has a maximum value 4/9 ~ 0.4444 at t = 2n/9 ~ 
0.6981. The naive solution has a maximum value 5/9 ~ 0.5556 at t — (10/9) 3,/2 ~ 1.1712. 
The renormalized solution has a maximum value 20/(27^) ~ 0.4277 at t = (20/27) 3 ^ 2 ~ 
0.6375. The renormalized solution gives accurate approximated solution and reproduces the 
contracting behavior of the universe qualitatively and quantitatively well. 

For T = 4/3, the energy density of the fluid is p\ oc a~ A and the effect of the fluid is the 
same as radiation. The naive and the renormalized solutions are 



a 



f 2/3 



i + (r 2/3 



a(t) 

l R (t)=t 2 / 3 ( c + 9ci r 2 / 3 ) 1/4 



A a 4 



-2/3 >, 



(15a) 
(15b) 



The renormalized solution behaves as a R oc t 1 ^ 2 for t ~ and a R oc t 2 / 3 for t ~ oo. Therefore 
the renormalized solution reproduces the behavior of the scale factor which represents the 
transition from the radiation dominant era to the matter dominant era. We cannot get such 
a behavior from the naive solution ( |15a| ). 
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III. ZERO MODE SOLUTION OF COSMOLOGICAL SECOND ORDER 

PERTURBATION 



We aim to treat the cosmological back reaction problem using perturbation approach. 
Let us assume that the metric is expanded as follows: 

9ab = 9ab + 9ab + 9ab + -- - ■ ( 16 ) 

(0) (1) 

g ab is the background FRW metric and represent the homogeneous and isotropic space, g ah 
is the metric of the first order (linear) perturbation. This metric represents the small linear 
fluctuation from the background space time. We can assume that the spatial average of the 
first order perturbation vanishes: 

( { 9ab) = 0, (17) 

(2) 

where (• • • ) means the spatial average with respect to the background FRW metric, g ab 
is the second order metric and contains non-linear effect caused by the first order linear 
perturbation. This non-linearity produces homogeneous and isotropic zero mode part of the 
second order metric. That is 

< <L) * 0. (18) 

As we want to interpret the zero mode part of the metric as the background FRW metric, 
we must redefine the background metric as follows: 

(0) (0) l/ (2) \ na\ 

9ab ► 9ab+ \ 9abl- ( 19 ) 

This is the back reaction caused by the non-linearities of the fluctuation and it changes 
the background metric. But in cosmological situations, the second order perturbation term 
grows in time and dominates the background metric, this simple prescription does not work 
well in the context of the perturbation expansion. Furthermore, the meaning of the gauge 
invariance is not obvious in the second order quantity. We cannot adopt Eq.(|T9"D as the 
definition of the background metric because the gauge transformation changes the definition 
of the background metric. We will resolve these problems in the next section by constructing 
the second order gauge invariant quantity and applying the renormalization group method. 

In this section, we obtain the solution of the second order perturbation. To extract the 
effect of the back reaction, it is not necessary to know the full form of the solution. We 
only need the homogeneous and isotropic zero mode of the second order perturbation. We 
consider the perturbation of the Einstein's equation with dust field in comoving synchronous 
gauge. The background is assumed to be a spatially flat FRW universe. The metric is 

ds 2 = -dt 2 + g i:j (t, x)dx i dx j . (20) 

The Einstein's equations are 

(3) i? + K 2 - K { jK\ = 2p, (21a) 

#4- -#,4 = 0, (21b) 

K t j +KK i j + ^ j = ^S i j . (21c) 
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where ^Rij is the Ricci tensor of the spatial metric, p is the energy density of the dust and 
the extrinsic curvature is defined by 

K i i = ^g jk g ki . (22) 
By eliminating p, the evolution equation can be written 

K\ + KK\ - 1 5h(K 2 - K l m K m {) + {^Rh - ± S j ^R^ = 0. (23) 

Using the conformally transformed metric 7^ = a~ 2 (t)gij, 

K\ = H8i i + l - 1 i l %, (24) 

where H = a/a and a is the scale factor of the background FRW universe. The evolution 
equation for the three metric becomes 

a V \ J 1 (25) 

= (2iJ + 3# 2 ) - - k + - ^ (A; 2 - fc' m k m t ), 

Z o 

where is the Ricci tensor of the conformally transformed metric 7^ and fc^ = 7^. We 
solve this equation perturbatively up to the second order: 

lij = Sij + hi j + l, r (26) 

where hy is the metric of the first order and 1^ is the metric of the second order, respectively. 

A. Background solution 

The background evolution equation is 

2H + 3H 2 = 0, (27) 

and the solution is 

a(t)=a t 2/3 , (28) 

where a is a constant of integration. 

B. First order solution 

The first order evolution equation is 

2 /W. 1 . (!) \ 

h?i + mh\ + -[ R\{h) - -5\R(h) \ = 0, (29) 
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(1) 

where Rij(h) is the linear part of the Ricci tensor: 

(1) 1 , i 

Rij(h) = 2 {-Kn - v h a + W + )• ( 30 ) 

The momentum constraint is 

K\i - \i = 0. (31) 
1. scalar mode 

For the scalar mode perturbation, the metric can be written using scalar functions as 
follows: 

hf = A(t,x)5 ij + B >ij (t,x). (32) 
By using the momentum constraint fl3~ip, we have 

A(t,x)=0. (33) 

Therefore 

2D 

A(t,x) = -*(x), (34) 

where $ is an arbitrary function of the spatial coordinate. The evolution equation becomes 

2D 

( V 2 BY+ 3H(V 2 B) ■ - —-.V 2 * = 0, (35) 

27 cr 



and the growing mode solution is 



2t 2 / 3 

fl=-5-¥. (36) 



The scalar mode solution is given by 

hf = ^*(x)dy + £L-^(x). (37) 
The density contrast is given by 



(i) 

V 2 ^ = V 2 ^ oc a(t). (38) 



J- = v * = 

(o) 9 a 2 H 2 a, 



p *o 
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2. tensor mode 



The evolution equation for the tensor mode is 

h^ W) + 3H h^ W) - 1 V 2 h^ W) = 0, (39) 

where h^f W ^ satisfies the transverse traceless condition: 

hi j = hi = 0. (40) 

The solution is given by 

fc ( 4i ) 

; A^t, k) = A tJ (k) r 1 ' 2 Z ±3/2 ^ t 1 ^ > 

where Z is a Bessel function and -Ay(fe) is a constant tensor satisfying the transverse traceless 
condition. 

C. Second order solution 

The evolution equation for the second order is 

..2 A 2 ) r . 1 • ( 2 ) r 
p. + 3H p. + _ I R LJ.__ 5 3. R L 

2 A 1 ) 1 (!)\ • • 1 • • 1 • • • 

- - - - 5 ki R + ^ fe h ki --h b?i + - ^ (h 2 - h\ h l k ) (42) 

er \ 4 / 2 8 



2 

^2 



i^L,. _ h jk Rk . _ _ 5 j. ( ^JVL _ ^ ^ 



(2) 

where Rfj is the linear part of the Ricci tensor with respect to the second order metric 

(2) 

and Rfj L is the quadratic part of the Ricci tensor with respect to the first order metric hij. 
Due to the homogeneity and the isotropy of the background metric, the zero mode part of 
the second order metric satisfies 

<^> = ^M0, (43) 
where the spatial average is defined as 

(li d ) = ^J<Pxl iS . (44) 

V is the volume of a sufficiently large compact domain and we assume periodic boundary 
conditions for perturbations. Therefore it is sufficient to consider only the trace part of the 



9 



second order metric. The trace part of the evolution equation ( f42|) with the spatial average 
becomes 

/ 5 . 1 • 2 /I ( 2 ) 1(1) 3 (!) \ \ 

(/)• + 3F <i>= (~ h» h* ~ \h* - 1 (\r nl - \h R + \h kl R k A }, (45) 



where 



(2) 1 

R NL — 



1 



h k i,i h kl ' 1 + h kl (y A h kl + h kl - 2 h kijl ' 



2 L2 

,»W /J, l. ^ / ukl ^ u,k 



(46) 



+ h^ 1 {h iKl - hu, k ) - [ h kl .i - - h> K ) (2 h k f - h kj 



1. second order zero mode solution : contribution of scalar mode 
For the first order scalar mode solution fl37|), Eq.fl4*5|) becomes 



i nn j.— 4/3 90 +—2/3 



where \l/fc is the Fourier component of ^f(x). The solution of this equation is 

i n /2/s +4/3 

(o = T — E e + — E fc4 **** + const - ( 48 ) 



This solution is consistent with the result of Tomita and Kasai p0|j21[| who obtained the 
complete form of the second order solution. 

2. second order zero mode solution : contribution of tensor mode 



For the first order tensor mode solution (41), Eq.(45) becomes 



(/)•■+ 3H (I)' = E (jj A k % k)A* kl (t, k) - ^ A k \t, k)AUt, k)j . (49) 

For the mode whose wavelength is smaller than the horizon scale k ^> aH, the first order 
solution is approximated to be WKB form: 

A tJ (t,k) « -A i:j {k)e ik ff , A\ = k l A tJ = 0. (50) 

(X 

For this solution Eq.fl4"9"|) becomes 

(/)•■+ 3H (/)• « --L E k 2 A^(k)A^(k), (51) 



4 a 4 
k 
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and the solution is given by 



q +-2/3 



(52) 



For the long wavelength mode whose wavelength is longer than the horizon scale k <C aH, 
the first order tensor mode solution becomes 



Aij(t, k) ~ Ajj(fc) 



and the Eq.(|49D becomes 



5\aH 



(53) 



(/)■■+ sh (iy « £ /c 2 ^(fe)A*,(fe). 



(54) 



The solution is 



(55) 







IV. GAUGE INVARIANCE AND RENORMALIZATION OF THE SECOND 

ORDER ZERO MODE SOLUTION 



The second order perturbation has the zero mode part which comes from the the non- 
linearity of the Einstein's equation, and this changes the evolution of the background FRW 
universe. We want to treat this back reaction effect by using the renormalization method. 
The background FRW solution has a constant of integration which corresponds to the free- 
dom of constant rescaling of the scale factor. By absorbing the second order zero mode part 
to this constant using the renormalization method, we can obtain the effective scale factor of 
the FRW universe which contains the effect of back reaction due to the inhomogeneity. But 
the existence of the gauge freedom prevents us from direct application of the renormalization 
method. If we change the gauge condition, the component of the second order zero mode 
solution may change. We do not know what component or combination of the zero mode 
solution should be renormalized. 

In this section, we first review the method of construction of the second order gauge 
invariant quantities using Abramo's argument ||-10|. Then we apply the renormalization 
group method to the second order gauge invariant quantity for the zero mode. Let us 
consider the infinitesimal coordinate transformation 



x 



(56) 



where £ M is assumed to be the first order small quantity. An arbitrary tensor q receives the 
gauge transformation 



q 



Q 



q = q- £$q 



(57) 
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where is the Lie derivative with respect to the vector £. Now we prepare the vector X 
which transforms under the coordinate transformation as follows: 

X» — >x'^ = X tM + ^ + ^[X, £}>* + ■■■ . (58) 

Then the quantity Q = e Cx q transforms as 

Q — >Q' = e Cx ' q' = e Cx, e~ c t q 

= e Cx q = Q. (59) 

Therefore the quantity Q is gauge invariant. If expand the tensor q perturbatively q = 

(0) (1) (2) 

q + q + q , the gauge invariant quantities in each order are 

(°) (0) 

Q = q, (60a) 

W (1) (0) 

Q = V + Cxq, (60b) 

/n\ / (2) \ , ir , 1 ir* (0) \ , ir (0) \ fan \ 

\Q) = \ q) + \£x q) + -{C x q) + {C Y q}. (60c) 

(0) 

Q is the background variable and does not receive gauge transformation. Y represents the 
vector which is needed to construct the gauge invariant combination under the second order 
coordinate transformation. We does not need explicit form of this vector. For the second 
order zero mode metric, the effect of Y is to produce the second order time coordinate 
transformation. We can set Y = if we do not mind the freedom of the second order 

(2) 

time coordinate transformation. The quantity (Q) is invariant under first order gauge 
transformation for Y = 0. 

Now we look for the explicit form of the gauge invariant quantity in a specific gauge. We 
consider the following form of the first order metric: 

(i) _ (-2<t> a B i \ 

9 ^ \aB ti a 2 (-2^5 i , + 2E, J )J " {0L> 
Under the coordinate transformation, the perturbation variables transform as follows 

0' = - i°, (62a) 
B' = B-a(\A' + -e, (62b) 



a 2 J a 

iP' = ^ + H^° } (62c) 

E , = E-\t, (62d) 
cr 

5x' = 5 X -eXo, (62e) 

where the scalar function \ represents the velocity potential for the dust field. From this, 
we can choose the following vector X which obeys the desired transformation (|58|): 

X = [-f-5 ij Ej\. (63) 
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Of course, the vector X is not unique. There are infinite possibility of the form of X which 
corresponds to the infinite possible candidate of the gauge invariant combination of the 
perturbation variables. We choose this form because we want to impose the comoving gauge 
condition 5\ = for the perturbation variables. The first order gauge invariant quantities 
constructed from X are 



(i) 
Qoo 

(i) 

Qoi 
(1) 

(i) 

Qy 



VXo 

R 2 TP I 6X 

an — a E H 

Xo 

H 

Xo 



5 X + X°xo = 0. 



-2 a 2 5 i:j ( ip + ^ 6x 



(64a) 
(64b) 
(64c) 
(64d) 



For comoving synchronous gauge <fi = B = 5\ = 0, there is no freedom of the coordinate 
transformation and the each component of the metric perturbation can be written using the 
combination of these gauge invariant variables. 

A. gauge invariant zero mode quantity: contribution of scalar mode 



Using the first order scalar mode solution (|37|), the zero mode gauge invariant quantities 
in comoving synchronous gauge are 



(2) (2) (2) 

(Qoo) = (Qo,) = (Q x )=o, 

Qij + {Qij) = a 5ij + {g ij) + {Lx 9 y > + - {£ x 9 ij) 



= a 2 5 



1 + ~ {ip, k E, k ) - i (E M E M ) 



L / (2) \ 
+ ( 9 ij) 



= a 



, 10(2/3 2, , 

fc 



Sij. 



(65) 



These quantities are invariant under the first order gauge transformation. To interpret the 

(0) (2) 

geometrical meaning of the gauge invariant quantity Qij + {Qij) } we consider the determinant 



of the three metric: 



(Vg) 



1 + E,) - L - (E# E ttj ) + f + \ (f) 



1 (0) (2) 

-^(Q^ + (Q^)) + a 2 (^ 2 )+const.) 



(2) 



2 

3/2 



'1 .. (0) (2) 

3*"(Q<, + <Q y » 



3/2 



(66) 
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(0) (2) 

Therefore, the trace of the gauge invariant variable + (Qy) corresponds to the physical 

(0) (2) 

volume element of the universe and we can interpret Q + (Q) as square of the scale factor 
of the FRW universe: 

ds 2 = -dt 2 + d 2 (t) dx 2 , 

; ^(M)) 1 ^ ^-gp*,,;). <67) 

At this stage, we renormalize the expression for the effective scale factor a(t) to improve its 
behavior for large t. The background FRW solution has the constant of integration ao- We 
regard the second order term in a(t) as a secular term and apply the renormalization group 
method. The renormalization group equation for ao becomes 

dt 2/3 27 a ^ k fe ' 1 ' 

k 



and the solution is 

I 2 

/ ° \ 
- \^k 2 tyuW I 

27 



a (t)= (c-^^A; 2 ^j , (69) 

where c is a constant of integration. The renormalized expression of the line element (^) is 

ds 2 = -dt 2 + (a R {t)) 2 dx 2 , 

Comparing with the analysis of the section II, this solution is the same as the FRW equation 
with positive spatial curvature. We conclude that the effect of the inhomogeneity due to the 
scalar mode fluctuation is equivalent to positive spatial curvature. This result is the same 
as Russ et al. who derived the spatially averaged Einstein's equation for a comoving 
volume and evaluate its behavior using the solution of the second order perturbation of the 
Einstein's equation. 



B. gauge invariant zero mode quantity: contribution of tensor mode 

For the first order tensor mode perturbation, the second order solution (Zy) is invariant 
under the first order gauge transformation because the first order tensor mode is already 
gauge invariant. The gauge invariant zero mode quantity is given by 

(0) (2) 

Qv + (Q tJ ) = a 2 (5 t3 + (h 3 )). (71) 
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Using the second order zero mode solution, we have 

3r 2 / 3 



(0) (2) 

Qa + (Qa) = { 



a 2 I 1 



a 2 I 1 



8 at k 



J2k 2 A ij (k)A* i:j (k) 5ij for k » aH 



21^3 

80 ~aj 



k 2 A ij (k)A* i:j (k) 5 lj for k<^aH 



(72) 



The renormalized scale factor becomes 



a R (t) = { 



1/4 



1/2 



for ^> aH 
for <C aH 



(73) 



Comparing with the result of section II, the effect of the back reaction of the short wavelength 
tensor mode is same as radiation fluid, and the long wavelength tensor mode is same as 
positive spatial curvature. This result is consistent with the analysis of Abramo [11] using 
the effective energy momentum tensor. 



C. gauge invariant zero mode solution in the longitudinal gauge 

To compare the result in the comoving gauge with other gauge condition, we use the 
longitudinal gauge E = B = here. We present only the result of the calculation. The first 
order solution for the scalar mode is 

9>*-\ -a 2 f^(x)5 lJ ) ' U4j 
and the second order zero mode solution is 

In this gauge, each components of the first order metric are gauge invariant quantities and 
we can set X = to construct the second order gauge invariant quantities. Therefore the 
gauge invariant zero mode solution is given by very simple form: 

(0) (2) (0) (2) 

Qij + Qij = 9ij + (9ij)- (76) 
This means the metric components of the line element 



k 2 * h *lj dt 2 + a 2 [l~^ 2 t 2/3 £ k 2 <M$J dx 2 (77) 
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are gauge invariant under the first order gauge transformation. As we have the freedom of the 
second order coordinate transformation, which is the second order coordinate transformation 
of time, we can transform this metric to the FRW form. We change the time coordinate as 

t = T- ! dt^-e'^^^l- (78) 



Then the metric becomes 



r/s 2 

""'long 



-dT 2 + T 4 / 3 ^a 2 - H T 2 / 3 ^ k " *fc*fcj dx '- (79) 



At this stage, we use the renormalization method and the second order perturbation term 
is absorbed to the constant a of the zeroth order. The renormalization group equation is 

^=4i5> 2 ****< < 80) 

k 

and the solution is 

/ \ 1/2 

a (T)=(c-^T^k^ k n) ■ (81) 

The renormalized metric is given by 

ds 2 ong = -dT 2 + (a R {T)) 2 dx 2 , 

. a R( T ) =T 2/3 a (T). (82) 

This metric give the same result in the comoving gauge. Therefore we have confirmed the 
gauge invariance of the back reaction effect. 



V. SUMMARY AND DISCUSSION 

We have investigated the cosmological back reaction problem using the second order 
perturbation of the Einstein's equation. To describe the back reaction effect due to the 
inhomogeneity in the framework of the cosmological perturbation approach, we used the 
renormalization group method. The second order zero mode solution which appeared by the 
non-linear effect is regarded as a secular term of the perturbative expansion, we renormalized 
a constant of integration contained in the background solution and absorbed the secular term 
to this constant. The renormalized constant obeys the renormalization group equation which 
determines the dynamics of the effective background universe. 

Owing to the gauge freedom of the Einstein's equation, it is not obvious that what com- 
bination of the second order variables should be renormalized. Constants contained in the 
background solution change if we use a N(t) ^ 1 lapse function for the background metric. 
The gauge independent background constants can be obtained by using the Hamilton- Jacobi 
method p2] and we can recognize the constant a , which is the freedom of rescaling of the 



16 



scale factor, has the gauge invariant meaning. For perturbation, using the second order 
gauge invariant variable, we renormalized the secular term to the background constant in 
a gauge invariant way. In comoving synchronous gauge, the gauge invariant combination 
corresponds to a comoving volume of the universe. This means that the requirement of the 
gauge invariance for inhomogeneous spacetime uniquely picks out the effective background 
FRW universe. In previous works, the effective scale factor was introduced through the 
uniform Hubble slicing condition |5],|7| and the comoving volume was introduced from the 
first ||. Anyway gauge invariance of the back reaction was not manifest. 

For the dust dominated universe, we derived the renormalization group equation and its 
solution is obtained. The scalar mode and the long wavelength tensor mode inhomogeneity 
work as positive spatial curvature, and the short wavelength tensor mode works as radiation 
fluid. These results are consistent with previous works PJ7|JTl[|. We also confirmed that 
these results are gauge independent by comparing the calculation using the comoving and 
the longitudinal gauge. Besides determining the dynamics of the effective background FRW 
universe, the renormalization group method gives the evolution of the fluctuation which 
includes the back reaction effect. As the background evolution is affected by the inhomo- 
geneity, the evolution of the fluctuation is also modified. For the first order perturbation 
solution of the density contrast fl58p, we can obtain the renormalized evolution by replac- 
ing do contained in the solution with the renormalized value ao(t). For the scalar mode 
perturbation, we have 

W £2/3 £2/3 

This indicates that the back reaction effect of the inhomogeneity enhances the growth rate 
of the density contrast. 

Including the scalar field is an interesting task to investigate the back reaction problem in 
the inflationary cosmology. This was partly done by Abramo et al. P^TH|] using the effective 



energy momentum tensor, but they do not derived explicit form of the evolution equation 
of the effective scale factor and the scalar field. We will treat this subject in a separate 
publication. 
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